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Breakdown of Fermi liquid behavior near the hot spots in a two-dimensional model: A 

two-loop renormalization group analysis 
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Motivated by a recent experimental observation of a nodal liquid on both single crystals and 
thin films of Bi2Sr2CaCu20g+i by Chatterjee et al. [Nature Physics 6, 99 (2010)], we perform a 
field-theoretical renormalization group (RG) analysis of a two-dimensional model consisting of eight 
points located near the "hot spots" on the Fermi surface which are directly connected by spin den- 
sity wave ordering wave vector. We derive RG equations up to two-loop order describing the flow 
of renormalized couplings, quasiparticle weight, several order-parameter response functions, and 
uniform spin and charge susceptibilities of the model. We find that while all the order-parameter 
susceptibilities investigated here become non-divergent at two loops, the quasiparticle weight van- 
ishes in the low-energy limit, indicating a breakdown of Fermi liquid behavior at this RG level. 
Moreover, both uniform spin and charge susceptibilities become suppressed in the scaling limit 
which indicates gap openings in both spin and charge excitation spectra of the model. We argue 
that all these results consistently suggest that the low-energy state of this model may indeed satisfy 
some key conditions to be characterized as an insulating spin liquid, giving further support to a 
previous analysis based on a one-loop RG calculation for this model available in the literature. 

PACS numbers: 74.20. Mn, 74.20.-z, 71.10.Hf 



I. INTRODUCTION 

The nature of the pseudogap state which manifests it- 
self in underdoped cuprates continues to generate both 
interest and controversy in the field of high-Tc supercon- 
ductivity. This is because the proper understanding of 
this phase turns out to be crucial for comprehending the 
underlying mechanism of superconductivity displayed by 
these materialai^. In recent years, angle-resolved pho- 
toemission spectroscopy (ARPES) shed new light on this 
problem. Experiments^ performed at very low temper- 
atures by Chatterjee et al. on both single crystals and 
thin films of Bi2Sr2CaCu208+5 gave important evidence 
pointing to the existence of an electronic state - located 
in between the antiferromagnetic (AF) and d-wave su- 
perconducting (SC) phases - whose excitation spectrum 
becomes zero only at the so-called nodal points (i.e. in 
the direction along the line connecting F = (0,0) and 
M = (tt, 7r)-points in momentum space). This state 
has thus been called a nodal liquid^"—. This electronic 
state is insulating but with a d-wave gap structure and 
exhibits non-Fermi liquid behavior around the so-called 
"hot spots" (i.e. points in momentum space where the 
AF zone boundary intersects the underlying Fermi sur- 
face of the system) . 

A minimal electronic model which potentially contains 
the rich phenomenology of the phase diagram displayed 
by the cuprates is the two-dimensional (2D) Hubbard 
model on a square lattice. This conjecture, originally put 
forward by Anderson^ in the mid 80 's, sparked among 
other things a profound reexamination of the founda- 
tions of Landau Fermi- liquid theory and, particularly, the 
precise conditions in which this standard theory breaks 
down in some 2D strongly correlated models. Indeed, 



soon after the experimental discovery of the cuprate su- 
perconductors, Anderson proposed that the ground state 
of the 2D Hubbard model slightly away from half- filling is 
given by an insulating spin liquid (ISL) which has no bro- 
ken symmetries down to zero temperature due to strong 
quantum fluctuations and also exhibits no well-defined 
quasiparticle excitations at low energiea^i^. This ques- 
tion, however, remains unresolved to this date. Despite 
that, these ISL states are very interesting from a theoret- 
ical viewpoint due to the fact that, once they are lightly 
doped with holes, they indeed give rise to a singlet d-wave 
superconducting phase at mean-field level^ in agreement 
with experimental observation. 

On the other hand, for a long time, there has not been 
a rigorous proof of the stability of this ISL state in re- 
alistic microscopic models. Part of this imbroglio owes 
to the fact that, at the present time, there is no gen- 
eral analytical or numerical technique that allows one 
to solve in a exact (or nearly exact) way such strongly 
correlated models in 2D. In the last years, important ad- 
vances have been made with improvements in both ana- 
lytical methods (see, e.g., Refs?^"— ) and also numerical 
technique&i^^— which, supplemented by groundbreaking 
experimental resultsi^Ti^i, provided strong evidence for 
the existence of these new quantum states in some real- 
istic physical situations. 

From a weak-to-moderate coupling perspective, renor- 
malization group (RG) methods remain one of the most 
powerful tools to attack these problems in view of its 
unbiased nature^^^^'*. They have been applied to many 
electronic models ranging from simpler 2D models^^^— 
which focus on the dominant role of scattering processes 
involving Fermi surface regions of the model to the fully 
2D Hubbard model defined on different lattice types. Re- 



garding the 2D Hubbard model on a square lattice for 
weak-to-moderate couplings, for instance, some of these 
RG works successfully reproduced both at one-loop^"— 
and two-loop orders'^'^''^- an AF phase in the model near 
half-filling, the onset of a d-wave singlet SC phase away 
from half-filling, and an additional electronic phase rem- 
iniscent of the pseudogap state^i^i^ which interpolates 
between these two phases. This agrees qualitatively with 
the physics displayed by the cuprate superconductors and 
gives further support to the point of view that this 2D 
model might indeed capture several important features 
of these strongly-correlated materials. 

Motivated in part by those theoretical findings, we 
present in this paper a two-loop field-theoretical RG cal- 
culation of the renormalized couplings, the quasiparti- 
cle weight, several order-parameter response functions, 
and the uniform spin and charge susceptibilities of a 
2D fermionic model consisting of eight points located 
near the "hot spots" on the Fermi surface which are 
directly connected by spin density wave (SDW) order- 
ing wave vector. A variation of this model was re- 
cently investigated in the literature using RG methods 
by Abanov and Chubuko\«2aiia, Metlitski and Sachdev^, 
and Efetov et als^, in which the corresponding high- 
energy fermions are partially integrated out in the sys- 
tem such that they arrive at a low-energy effective the- 
ory involving the bosonic SDW order parameter coupled 
to the fermionic excitations near the "hot spots". They 
found interesting renormalizations of several parameters 
of the model and also a clear breakdown of Fermi liq- 
uid behavior near these Fermi points. By contrast, our 
approach in the present paper follows closely the con- 
ventional RG strategySS, in which high-energy degrees 
of freedom are successively integrated out to derive an 
effective description in terms of renormalized low-energy 
fermion. This approach can be found, for instance, in the 
analysis performed by Furukawa and Rice^^ who first dis- 
cussed this model within a purely fermionic RG scheme 
up to one-loop order. As a consequence, these authors 
found, for weak Hubbard-like initial conditions, a RG 
flow of the renormalized interactions towards strong cou- 
pling at low energies, which seems to suggest some limi- 
tations of the one-loop RG scheme in order to access the 
infrared regime of this system. In our view, this fact also 
points to the importance of higher-order quantum correc- 
tions in the full description of the low-energy dynamics of 
this 2D model. In this respect, our present work takes se- 
riously this observation and, for this reason, it represents 
a step forward in this direction. 

This paper is organized as follows. In Sect. II, the 2D 
fermionic model consisting of eight points located near 
the "hot spots" on the Fermi surface that we wish to 
study is introduced. Next (in Sect. Ill), we explain the 
field-theoretical RG methodology and show how to im- 
plement this method up to two loops in order to describe 
this model. In this part, we will choose to present this 
methodology in a somewhat "condensed" way since the 
field-theoretical RG scheme up to two-loop order was al- 
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FIG. 1: The FS of the 2D fermionic model consisting of eight 
points located near the "hot spots" on the Fermi surface which 
are directly connected by spin density wave (SDW) ordering 
wave vector. 



ready described at length and in full detail in the con- 
text of another 2D fermionic model elsewhere by one of 
us*'^. The RG flow equations are derived analytically and 
solved numerically in Sect. IV and then we proceed to 
discuss our main results. Lastly, Sect. V is devoted to 
our conclusions. 



II. MODEL 

We begin our analysis with a general two-dimensional 
energy dispersion given, for instance, by ^k — 
—2t{cos{kx) + cos(fcy)) -t- At' cos{kx) cos{ky) — fi with t 
being the nearest neighbor hopping, t' the next-nearest 
neighbor hopping amplitude and fi the chemical poten- 
tial. For the cuprates, the appropriate choice of parame- 
ters is t' — 0.3i, which results in the curved Fermi surface 
(FS) shown in Fig. 1. This FS intersects the antiferro- 
magnetic Brillouin zone at eight points (i.e., the "hot 
spots"). These points are connected through Umklapp 
processes which leads us to examine the RG equations 
for all the coupling constants near these points. 

The eight points around the "hot spots" are displayed 
in Fig. 1. If we rotate the momentum axes {kx,ky) by 
45° degrees, we can define the new axes {k^^,k±), where 
the momenta fc|| and k± refer to momentum parallel and 
normal to the FS, respectively. Besides, since we will be 
interested only in the universal quantities of this model, 
we shall linearize the energy dispersion around the FS as 
Ck ~ |vf(^j|)|(|^_l| — kp) with the Fermi velocity vector 
given by VFik\\) = \'^{k^^,k^)^i<:\k^=kir\r>., where kp is the 
Fermi momentum at the hot spots and n is a unit vector 
perpendicular to the FS. In terms of the parameters of 
the two-dimensional energy dispersion, the absolute value 
of the Fermi velocity is given by 
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where A = IGt^ cos2(fc||/y2) + 16t'[l + cos2(/cj|/\/2)][/i + 

4t' sin^(fc||/\/2)]- Both the momenta paraUel to the FS 
and perpendicular to the FS are restricted to the interval 
[—kc,kc\, where k^ essentially determines the ultraviolet 
(UV) momentum cutoff in our theory. This implies also 
an energy cutoff which is given by Aq = 2vFkc which we 
choose to be equal to the full bandwidth of problem, i.e. 
Ao = 8t. 

If we use a coherent-state functional integral repre- 
sentation of the resulting Hamiltonian, the model at 
T = and constant chemical potential fi — Ep 
becomes described by the partition function Z — 
J Vlip, ip] exp(i J^ dtL[tp, tp]) with the Lagrangian L — 
Lq + Lint given by 



k,CT 

-E E 5.,B?.(k4)?<,,(k3)V.'(k2)^.(kl) (3) 

i ki,k2,k3 

cr,(T 

where k4 = ki + k2 — ka and the volume V has been 
set equal to unity. The Grassmann fields V'o-(k) and 
ipaO^) are associated, respectively, to the creation and 
annihilation operators of excitations lying in the vicin- 
ity of the hot spots with momentum k and spin pro- 
jection a. The index i runs over all possible inter- 
action processes of the model that produces logarith- 
mic divergences within perturbation theory, i.e., i = 
1, 2, 3, Ic, 2c, Ix, 2x, Is, Ir, 3x, 3p, 3t (for details of all the 
couplings taken into account, see Fig. 2). In this way, 
in order to keep a close connection with other RG works 
in the literature, we are following a "g-ology" notation, 
adapted of course to our 2D problem at hand. The La- 
grangian of Eq. ^ therefore defines our bare quantum 
field theory model which is regularized in the UV by the 
cutoff kr mentioned above. 



III. TWO-LOOP RG METHODOLOGY 

The methodology of our RG scheme follows closely the 
standard field-theoretical approach^, which was also ex- 
plained in full detail in a previous paper published by one 
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FIG. 2: The relevant couplings of the field theory model. 
Here we use a "g-ology" notation adapted to our 2D problem 
at hand. 



of ue"*"^ . If one applies a naive perturbation theory for the 
present model, divergences and other non-analyticities 
emerge in the low-energy limit at the calculation of sev- 
eral important physical quantities of the model such as 
two-particle effective interactions, the self-energy, and 
several susceptibilities. This result implies that the con- 
ventional perturbation theory setup is not appropriately 
formulated for this case. We circumvent this problem 
by defining renormalized quantities and their appropriate 
counterterms in the theory which will effectively regular- 
ize the model at a fioating RG scale parameter A. As a re- 
sult, the new renormalized perturbation theory becomes 
well-defined in the low-energy limit A — >■ 0. This regu- 
larization procedure has to be implemented consistently 
order by order in perturbation theory. In view of the fact 
that this program is successfully accomplished here, the 
present field theory model turns out to be renormalizable 
up to the order within perturbation theory calculated in 
the present work. 

At two-loop RG level, we must perform the following 
substitutions for the fermionic fields in the Lagrangian 
(Eq. (131)) of the model 



V^,(k) ^ ZyVfl,(k), 



^.(k) 



A/2' 



^fl<T(k), 



(4) 



where Za is the RG flowing quasiparticle weight which 
is naturally related in the limit of A ^ to the con- 
ventional many-body deflnition of the quasiparticle peak 
Z = {1- 5ReS(w, k = kF)/9a;|„=o)~\ with E(w, k) be- 
ing the self-energy and kp is the Fermi vector. The cor- 
responding Feynman diagrams up to two-loop order are 
shown in Fig. [31 The one-loop diagram (i.e. the Hartree 
term) is generally independent of the external frequency. 
As a result, this contribution does not renormalize the 
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FIG. 3: Feynman diagrams for the self-energy of the model up 
to two-loop order. The one-loop contribution is the Hartree 
term and the two-loop contribution is the so-called sunset 
diagram. 



quasiparticle weight in the present case. In fact, it only 
generates a constant shift in the chemical potential of 
the model that must be appropriately subtracted by a 
counterterm such that the density of particles in the sys- 
tem remains always fixed^^. By contrast, the two- loop 
contribution (i.e. the sunset diagrams) is the first contri- 
bution to the self-energy that produces a nonanalyticity 
as a function of the external frequency oj. For this reason, 
this term alone will be responsible for the renormaliza- 
tion of the quasiparticle weight up to this order (for more 
details on this point, see, e.g., Refs4^i^ in the context of 
different fermionic 2D models). 

The so-called anomalous dimension is conventionally 
defined by 77 = Adln^A/dA (see, e.g., Ref4i). From this 
expression, we obtain using the standard RG condition 
for the inverse of the renormalized single-particle Green's 

.(2) 



function r]^, (fco ^ A,kj_ = kp) = A that 
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Here we would like to stress the fact that if we were to 
take the ID limit of the above equation, we would re- 
produce exactly the well-known result for the anomalous 
dimension for this case at two-loop order as was first cal- 
culated long ago in Ref."*^. 

In addition, it can also be conventionally defined^ that 
the bare and renormalized coupling constants are related 




FIG. 4: Feynman diagrams for vertex corrections up to two- 
loop order. The last diagram represents the counterterm. 



to each other at two loops by 
9i.. 



N-\0)ZX' 



g..n + ^9^' + ^aT'' 



(6) 



where i = 1, 2,3, Ic, 2c, la;, 2a;, Is, Ir, 3a;, 3p, 3t and the 
term A^(0) = kc/{TT'^vp) is the density of states of the 
model at the "hot spots" points. The renormalized quan- 
tities - labeled by the subscript R - generally depend 
on the RG scale A. In contrast, all the bare quantities 



Hoop 



and 



will be denoted by the label B. The terms Ag^j^ 
^9iR°^^ represent the counterterms necessary to regu- 
larize at one-loop and two-loop orders, respectively, the 
two-particle irreducible vertex corrections in each of the 
corresponding scattering channels. 

We can now adjust these counterterms such that all di- 
vergences are exactly canceled in our series expansion for 
the couplings up to two loops, i.e. we choose the following 
standard RG condition for the two-particle irreducible 
vertex functions: r^^(fcio = A/2,fc2o = A/2, fcao = 
3A/2,fc4o = —A/2) = —igi^]i{A), where kio are the en- 
ergy components. However, the price we pay for this is 
the appearance of a new scale with all physical quan- 
tities now depending on this scale A. By contrast, the 
original model has no information about this quantity, 
i.e. the bare parameters do not depend on A. This leads 
us to the renormalization group conditions for the bare 
couplings of the model, i.e. KdgiB/dK = 0. Therefore, 
using Eq. ([B]), we finally obtain 
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Agf^r' 



(7) 



The initial conditions for this system of differential equa- 
tions are naturally given by the microscopic interactions. 
A representation of the Feynman diagrams correspond- 
ing to the vertex corrections up to two-loop order are 
displayed in Fig. ID 

In order to investigate what are the enhanced correla- 
tions in the low-energy limit of the model, it is important 
to calculate the corresponding susceptibilities by intro- 
ducing an infinitesimal external field in the appropriate 
channel and evaluating its linear response. Therefore, we 
must add to the Lagrangian that describes the present 
model the following term 



E'^B,SC^a(k)^^(-k) 



k,Q,0 



E rs' 



DW 



^„(k + Q)VMk) + i^.c., 



(8) 



k,a,0 



where Tb sc ^'^^ ^b dw ^^'^ ^^'^ bare response vertices 
for superconducting (SC) and density-wave (DW) orders, 
respectively, and Q = (7r,7r). This added term will gen- 
erate new Feynman diagrams - the three-legged vertices 
displayed in Fig. [5] - which will also produce new loga- 
rithmic singularities in the low-energy limit of our field 
theory model (see also Refsi^i^ in the context of differ- 
ent fermionic 2D models). Therefore, we must regularize 
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FIG. 5: Feynman diagrams for both SC and DW response 
functions of the model. The diagrams with crosses represent 
the counterterms. 



these divergences by defining the renormahzed response 
vertices and the corresponding counterterms as follows 
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As a result, by invoking the RG condition for bare quan- 
tities of the model K{dTg gQ,jjy^r)/d^) = 0; we obtain 
the RG equations for the SC response vertices 
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where we defined the new vertices T^ 

Vi!scih^^^) ^^'^ tI^SC = '^,sc(^lh^^-L) depending 
on their location at the Fermi points. As for the DW 
response vertices we get the following RG equations 
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where j = 1,2 and the subscripts SSC and TSC corre- 
spond to singlet and triplet superconductivity, whereas 
SDW and CDW stand for charge and spin-density 
waves, respectively. Hence, using the above relations one 
can derive in a straightforward way the RG flow equa- 
tions for each response vertex associated with a poten- 
tial instability of the normal state toward a given ordered 
(i.e. symmetry-broken) phase. In order to determine the 
symmetry of the order parameter we must further sym- 
metrize the response vertices w.r.t. the index i. Thus 
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Here the labels s — wave, d — wave and p — wave refer 
to the symmetries of the corresponding SC order param- 
eters, and CF and SF stand for charge and spin flux 
phases, respectively. 

We now move on to discuss the uniform (q — >■ 0) linear 
response functions of the model. At one-loop level, we 
find that there are no Feynman diagrams that produce 
logarithmic singularities in the low-energy limit. In this 
way, it becomes necessary to go up to two loops in order 
to obtain the appropriate RG equations for this quantity 
(see Fig. [6]). To do this, we must add new bare auxil- 
iary fields T§'"(q« 0) that couple, respectively, to the 
charge and spin densities in the Lagrangian of the model 
as follows 

L'ext - E [ri^(q)?t(k)^T(k) + ri^(q)?;(k)^;(k)] . 

k 

(12) 



where we also introduced the new vertices T^ j^^/ = 

TRjDwih ' ^-l) and T^^^j^^ = T^jjwih ' -^J-) according 
to their association to the hot spots. By symmetrizing 
and antisymmetrizing all the response vertices w.r.t. the 
spin indices, we obtain the following order parameters 



Thus, we must regularize these divergences by defining 
the renormalized uniform response functions and the cor- 
responding counterterms as foUows^i^ 



mcL) = ZX' [Tri<D + Ar«""(q)] , 



(13) 
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FIG. 6: Feynman diagrams for uniform (q — >■ 0) response 
functions of the model up to two loops. Since the correspond- 
ing one-loop diagrams for this case are non-divergent near the 
hot spots, they turn out not to be relevant in this calculation. 
The last diagram stands for the counterterm. 



where the Za factor originates, as discussed previously, 
from the redefinition of the fermionic fields at two-loop 
RG level displayed in Eq. @ and is naturally related 
to the self-energy feedback into the equations. In order 
to determine these counterterms, we must change the 
RG condition of the three-legged vertices to Fqq ' (go — 
0,VF{\q±\ — kF) = A) = — i7^"(A) since we will approach 
now the low-energy limit of the model via the momen- 
tum scale. Following this prescription, by establishing 
the RG condition for the bare quantities of the model 
A{dTg"/dA) = 0, we obtain the RG equations for the 
renormalizcd uniform response functions 



AS^ = „rr(.) 



dA 



By symmetrizing and antisymmetrizing these quantities 
w.r.t. spin indices, i.e. 



Tc = n' 
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tt 



r, 
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we readily obtain the corresponding fiow equations for 
both renormalized charge (Tc) f^nd spin (Ts) uniform 
response functions. 

Lastly, once we computed the response vertices asso- 
ciated with all the order parameters, we can proceed to 
calculate their corresponding renormalized susceptibili- 
ties. They are given by 



xm{A) = N{o) / d^r™(or^(o, 



(15) 



where m = SSC{s — wave), SSC{d — wave),TSC{p — 
wave),CDW,SDW,CF and SF, and the RG step / is 



defined by I — ln(Ao/A). As for the uniform suscepti- 
bilities, one obtains that x„(A) = iV(0)|7^i(A)p, where 
n = C and S. 



IV. TWO-LOOP RG FLOW EQUATIONS AND 
NUMERICAL SOLUTION 

In this section, we discuss the RG flow equations up to 
two-loop order associated with all the physical quantities 
calculated previously in this paper and show their nu- 
merical solution as one approaches the low-energy limit 
of the model. Concerning the renormalized couplings, we 
choose to write down for completeness all the RG equa- 
tions up to two-loop order explicitly in Appendix A. It is 
important to stress here that, up to one-loop order, our 
equations agree with those derived by Furukawa and Rice 
(Ref.^'^) who first discussed this model within a one- loop 
RG approximation in the literature. However, as we have 
explained previously in this paper, to go beyond the work 
in Ref.^* we need to investigate the effect of higher-order 
quantum fiuctuations in order to describe their impact 
on the low-energy dynamics of the system. 

To solve all the coupled differential RG equations nu- 
merically, we apply standard fourth-order Runge-Kutta 
method. As an initial condition for these equations, 
we take Hubbard-like repulsive parameters, i.e. g^ = 
[kc/{7r'^VF)] U for all couplings. For simplicity, we will 
always choose g\ = 0.5 in the present work. Using 
the experimental data from Ref/^, since the nodal liq- 
uid is stabilized for 4% of hole doping, this implies that 
vf ~ 1.89i at the "hot spots" and the Hubbard local 
interaction parameter U ~ 4.4i. We mention here, how- 
ever, that our results do not depend critically on this 
initial choice of the coupling constants and we emphasize 
that they are robust within an appreciable range of these 
physical parameters. 

In order to establish a direct comparison between the 
one-loop RG results and the two-loop RG data for this 
2D model, we first present a plot of the fiow of the cou- 
plings up to one-loop order only. As a result, the corre- 
sponding one-loop RG flow as a function of the step / is 
displayed in Fig. [T] In agreement with Refi^, we confirm 
that despite the fact that the couplings are initially taken 
to be weak, the one-loop RG fiow is always to a strong 
coupling regime in the low-energy limit for Hubbard-like 
initial interactions. Additionally, we observe that many 
renormalized couplings tend to diverge as the RG scale 
is lowered at the same critical scale Ic = ln(Ao/Ac), af- 
ter which point the RG approach up to one loop clearly 
breaks down. This fact seems to limit the validity of the 
one-loop RG calculation in order to study the infrared 
regime of this model and invites one to go beyond that 
approximation in order to access the low-energy limit of 
the system. 

Following this strategy, we now move on to our two- 
loop RG results. The RG fiow up to two-loop order for 
the couplings as a function of the step I is depicted in 
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FIG. 7: (Color online) One-loop RG flow for all the coupling FIG. 8: (Color online) Two-loop RG flow for all the coupling 
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Fig. [51 In this plot, we observe that instead of dis- 
playing a divergent behavior as shown in the one-loop 
RG flow, all renormalized couplings at two-loop RG level 
now clearly approach asymptotically infrared stable fixed 
points in the low-energy limit. Therefore, quite surpris- 
ingly, the two-loop RG approach turns out to yield a 
better controlled theory than the one-loop RG scheme 
for this model in the infrared regime and, for this reason, 
we can hope that it could describe (at least qualitatively) 
the correct trend of the low-energy dynamics of the sys- 
tem. Two increasingly relevant renormalized couplings 
in the RG flow turn out to be the Umklapp scattering 
g^x and the forward scattering 52 ■ As we will see shortly, 
they both tend to favor SDW antiferromagnetic ordering 
tendencies in the system. By contrast, upon inclusion of 
quantum fluctuation effects, the Cooper pair interaction 
processes 51c and 52c - that are intially taken to be re- 
pulsive - flow naturally to attractive infrared stable fixed 
points in the low-energy limit. These latter couplings in 
turn enhance SC pairing correlations which will of course 
tend to manifest itself as a possible competing order in 
the system. As can also be inferred from Fig. |S1 the 
remaining couplings flow asymptotically to zero in the 
infrared regime and, as a result, they become danger- 
ously irrelevant in the low-energy effective description of 
this 2D model. 

Next, we focus on the RG flow equation up to two- 
loop order for the quasiparticle weight Za of the model 
(see Eq. ([5])). By solving this equation numerically with 
initial condition Z\{1 = 0) = 1, we obtain the results de- 
picted in the plot of Fig. [5] as a function of the RG step I. 
We observe that this quantity becomes universally sup- 
pressed as a power-law Za ^ (A/Aq)'' with 77 « 2.01 in 
the low-energy limit, indicating a complete absence of 
coherent quasiparticles and a clear breakdown of Fermi 
liquid behavior near the "hot spots" for this regime. This 
conclusion is also in qualitative agreement with other 
works available in the literature^"—. Thus, in order to 
investigate the nature of this low-energy state, we must 



further characterize the model by examining what are the 
dominant fluctuations that drive the system to a nullified 
Z\ in the low-energy limit. 

For this reason, we now turn our attention to the be- 
havior of the various order-parameter susceptibilities as a 
function of the RG scale A of the model. Using Eqs. ^TU\\ 
and (jlip , we write down explicitly the corresponding RG 
flow equations for these quantities up to two-loop order 
as follows 



r|--«- = _(g2, + 5i, + 5i,+52x + 25i, + 2gi,) 



'ssc 



^s—wave 1^ ^s—wave 



'SSC 



'SSC ' 



1 



'T'ssc'"' = 2^-'^2c + 5ic-5i. -52. -251s + 2(71 

^ ^d—wave , ^d—wave 
^ 'ssc "•" '/ 'SSC ' 



1 



■p—wave 



'TSC — 2 ^^2^= ~ ^^'^ ^2^ ~ ■^1^-' TSC 

'^ '/ 'TSC ' 
TCDW = 2 (^51 + -^^ + ^^1^ + ^53p + 853* 



(16) 
(17) 
(18) 



93x- 92- g2x)TcDW +TlTcDW, (19) 



TsDW = ■:^i92+ g2x + 93 + 93x)TsDW 

+ vTsDw, 

%F = 2 (^51 + .92 - 2gij; + 52a; + 53 

- '293p- 93x)TcF +'nTcF, 

TsF = --^i93 + 93x- 92- g2x)TsF 

+ vTsF, 



(20) 
(21) 
(22) 



where we are using the shorthand notation 7i — 
A{dTi)/{dA). We solve these differential equations us- 
ing the same numerical procedure explained above and, 
by using Eq. (fT51) . we are able to follow the RG flow 
of the corresponding susceptibilities as a function of the 
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FIG. 9: (Color online) Two-loop RG flow for the quasiparticle 
weight Za of the model for the initial choice of g\ 
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RG step I. The results are plotted in Fig. [TOl In this 
figure, we obtain that even though many susceptibilities 
are renormalized due to interactions, they all remain non- 
divergent and become saturated at plateau values in the 
low-energy limit at two loops. In other words, the self- 
energy feedback and the higher-order vertex corrections 
at two-loop RG level are sufficient to suppress the diver- 
gence of these physical quantities displayed in one-loop 
RG calculations. Our present result implies that there 
should be no spontaneous symmetry breaking associated 
with any of the channels investigated in this work and 
therefore the system should not exhibit any long-range 
order at low energies. Indeed, despite the fact that the 
leading short-range correlations ~ by at least two orders 
of magnitude ~ for this regime are of SDW-type (i.e. an- 
tiferromagnetic spin correlations), spin flux phase fluc- 
tuations also appear as a subleading ordering tendency 
and d — wave singlet SC clearly comes in third place. 
This suggests that a possible coexistence of such short- 
range ordering tendencies could be the hallmark of the 
true nature of the corresponding low-energy state. 

Other physical quantities which are naturally of high 
importance to potentially identify the nature of the 
ground state of the model are the uniform (q — > 0) 
susceptibilities in both spin and charge sector. After a 
lengthy but straightforward calculation, we obtain that 
the resulting RG flow equations for the uniform response 
functions up to two-loop order are the following 



Ts 
fc 



1 



Ir 2 
2 [.93 



2 

9ix 



glp 



2 

9ic 



alx 



(gap- 93xf]Ts, 

{g3p-g3xf]Tc, (23) 



where % — A{dTi)/{dA). We emphasize here the fact 
that if we take, as an important check, the ID limit of 
the above two-loop RG equations for the uniform spin 
and charge response functions, we are able to reproduce 
established results available in the literature for these 
systems in various regimes (see, e.g., Ref*^). 



In the present 2D model, we once more follow the nu- 
merical procedure described before in order to solve the 
equations given by Eqs. (|23l) and Eq. (fTSJ) . Our result 
for the uniform charge susceptibility is shown in Fig. 1111 
We observe from this plot that, for repulsive Hubbard- 
like initial conditions, the uniform charge susceptibility 
(or simply charge compressibility) becomes strongly sup- 
pressed and renormalizes to zero in the low-energy limit. 
This implies that there should be a gap opening in the 
charge excitation spectrum in the system. This charge 
gap is produced by the Umklapp coupling parameter g^x, 
which, as we have seen before, becomes increasingly rele- 
vant in the low-energy effective theory. At this point, we 
draw attention to the fact that our result bears some re- 
semblance to the well-studied case of the two-chain repul- 
sive Hubbard ladder exactly at half-flUing'^"'^^, where it 
is by now well-established that the mechanism for charge 
gap generation in that system is driven by Umklapp inter- 
action in the scaling limit. Therefore, we may conclude 
from this that the nature of the low-energy state of the 
present 2D model is clearly insulating, which also agrees 
with our previous result that the quasiparticle weight Z\ 
is nuUifled in the low-energy limit. 

In Fig. [Til we also display the RG flow for the uniform 
spin susceptibility of the system for repulsive Hubbard- 
like initial conditions. As a result, we observe from this 
plot that even though this quantity is initially less sup- 
pressed in the RG flow than the charge compressibility, 
it also scales down eventually to zero in the low-energy 
limit. This indicates that there should be also a gap in 
the spin excitation spectrum of the model. The opening 
of a spin gap is produced by an interplay between the 
Umklapp coupling parameter g^^ (which also promotes 
SDW antiferromagnetic fluctuations and leads to charge 
gap formation) and the Cooper pair interaction process 
gic, which in turn renormalizes to an attractive infrared 
stable flxed point at low energies and enhances pairing 
correlations. This result strongly hints at short-range an- 




FIG. 10: (Color online) Two-loop RG flow for various order- 
parameter susceptibilities of the model (in units of A''(0)) for 
the initial choice of gl — 0.5. 







1.5 2 2.5 

/ (RG step) 

FIG. 11: (Color online) Two-loop RG flow for both charge (C) 
and spin (S) uniform susceptibilities of the model (in units of 
iV(0)) for the initial choice of gf = 0.5. 



tiferromagnetic spin correlations coexisting with singlet 
pair formation as the underlying physics of this corre- 
sponding low-energy state. 

Hence, it seems reasonable to suppose that our two- 
loop RG data points to the fact that the low-energy state 
of the present 2D model may indeed share some simi- 
larities to the d-Mott phase of the half-filled two-chain 
repulsive Hubbard ladder where the spin gap is well- 
understood to be related to a precursor effect of "pre- 
formed" Cooper pairs with approximate d — wave sym- 
metry in an otherwise insulating system^i^. This inter- 
pretation is in direct agreement with Ref.— , where the 
authors also concluded, based on a one-loop RG calcula- 
tion, that the low-energy state of this 2D model should 
be indeed given by an ISL state. However, we point out 
here that, within a RG computation up to one-loop or- 
der only, one cannot be sure about the spin properties 
of this state since the uniform spin susceptibility and the 
SDW spin susceptibility actually behave in a contradic- 
tory manner^^. Our present two- loop RG results resolve 
this contradiction in a most natural way. Since the SDW 
susceptibility at two-loop RG level becomes saturated at 
plateau values at low energies, it cannot be associated to 
any long-range spin order. This is also consistent with 
the fact that the low-energy state naturally displays a 
spin gap. Therefore, we believe that such an ISL state 
could be indeed produced by quantum disordering effects 
induced by a combination of coexisting short-range cor- 
relations included in the present 2D model, and the uni- 
versal properties of this state is captured qualitatively 
by our two-loop RG scheme within a weak-to-moderate 
coupling regime. 



ordering wave vector. By following the field-theoretical 
RG strategy, we have derived the corresponding flow 
equations for the couplings, the quasiparticle weight, the 
order-parameter susceptibilities and the uniform spin and 
charge susceptibilities of the model up to two-loop order. 
As a result, we have found for physically relevant choices 
of parameters that while SDW antiferromagnetism, spin 
flux phase, and d-wave superconductivity manifest them- 
selves as coexisting short-range ordering tendencies, the 
quasiparticle weight Z/^ renormalizes to zero in the low- 
energy limit, indicating the breakdown of the Fermi liq- 
uid behavior for this regime. Moreover, both uniform 
spin and charge susceptibilities scale down to zero in the 
same limit, which points to the fact that gap openings 
should should take place in both spin and charge excita- 
tion spectra of the model. By comparing those results 
with other well-studied systems such as the half-flUed 
two-chain repulsive Hubbard ladder, we have given ar- 
guments that our two-loop RG data are consistent with 
an interpretation that the low-energy state of the present 
2D model may have an insulating spin liquid-character 
displaying no long-range order at zero temperature. 

On the other hand, we would like to emphasize here 
that, despite the fact that the present 2D model might 
have some analogies with the well-known physics of the 
pseudogap phase of the cuprate superconductors, this 
model must be of course extended in order to describe the 
full Fermi surface (and, consequently, its process of be- 
coming partly truncated in momentum space) displayed 
by these compounds. For this reason, we believe that 
the present two-loop RG result could potentially give 
some insights to this problem starting from a weak-to- 
moderate coupling perspective. Lastly, we point out that 
it would be naturally very interesting to perform a com- 
plete two-loop RG analysis to discuss the fully 2D t — t' 
Hubbard model on a square lattice with a choice of pa- 
rameters that precisely match the Fermi surface observed 
experimentally in these materials. This will introduce 
some expected complications such as the substantial in- 
crease of the renormalized couplings as they become in 
this case functions of three independent momenta of the 
interacting particles at low energies, but this can be han- 
dled numerically using standard techniques. Therefore, 
we plan to perform such a RG investigation of the 2D 
t — t' Hubbard model in an upcoming publication^. 
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V. CONCLUSIONS 

We have performed a RG calculation of a 2D model 
consisting of eight points located near the "hot spots" 
on a Fermi surface which are directly connected by SDW 



Appendix A 

In this appendix, we show explicitly the RG flow equa- 
tions up to two loops for the model. They read 
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A ^9l 22^22 ^/ 2 \ 

A-^ = 91+ gix + 4g3t + 53p - 9ixg2x - gspgsx + ■^(gixg2x - g2x - gspgsxjgic 

+ 2^gic + gi+gix + g2x- gixg2x- g3pg3x + glp + glJgi, (Al) 

^IX "" 2 *^^^ ^ ^23; - 53 - .9L) + ^ (ff? + sicffL + gi9ic) + 1 (252 - gi) gl 

+ ^ [(2.g2c - fflc) (Ssp + SL) - '2g2cg3pgZx + 2^2 (g^p + SL - SSpffSa;)] 

+ ^{g2-g2c){glx + glx-gixg2x), (A2) 

^"dA ^ '■■^^ ~ ^^^-"^^ ^ -93* - gixigsx - '^gsp) - g2xig3p - 33^) + -[(.gi - 252)^ 

+ (51c - 252c)^ + 2gl^ + 2gl^ - 2gixg2x - ^gapga^ + 2g^^ + 2.g^^ + gl]g3, (A3) 

A-^ = glc + 9ixg2x + gls + glr + glx - gspgsx + -^ {gix92x - gL - gspgax) gi 

7^ {glc + gl + glx + glx - gixg2x + glp + glx - gspgsx) gic, (A4) 



2 
A-7T- = 1^ [gic + gix + g2x + gis + gn + ^gir - g-ip) + t \g\g\x + .9ic5i + g\. 



«<72c l,'2,2,2,2,2,r,2 2^,l/' 2 

7ic + g\x + g2x + 51. + .gi; + ^g\r 



dh 2 



-T ^(52c - g2){glx + SL - 5l:r52a;) + ^ (252c - 5lc) ffs 

+ 4 [(2.92 - ffi) (Ssp + sL) - '^g2g3pg3x + 252c (sfp + sL - Sspffsa;)] , (A5) 

A ^9ix , , „ 1 / , „ alp gix\ 

A-jT- = .gic.g2x + g2cgix + Agi^gir + ^ I .gig2c + gic.g2 - ^.g2cg2 — ^ ^ 1 51a: 

+ 27751,, (A6) 

A-T7^ = 5ic.gi2; + g2cg2x + "^gisgir + 2 ( gicgigix - 2gicgig2x + gicg2g2x 

+ gig2cg2x - 2g2cg2g2x - gixg3pg3x + g2xg3pg3x - -^g2xiglx + glp) j + '^vg2x, (A7) 

A-^ = {gic + g2c)gis + {gix + g2x)gir + '2vgis, (A8) 

A-TT" = {gic + g2c)gir + -^{gix + g2x)gis + '^v9ir, (A9) 

A-T^ = (251 - g2c)g3p + gixg3 + ^glt ~ g2g3p ~ 52^33 - 51532; + -^ f 252c5253p 

2 \ 

+ 5L33a: - 5l52c53p - 5lc5253p - gixg2xg3x - 5lc5l532; ^^^ j + 27753p, (AlO) 



A-TT- = (251c - 52c)53:e - 5lc53p - 52a;53 " 5253x + ^ ( 252c5253x + 52xff3p " 5l52c53a; 

2 \ 

- 5lc5253a; " 5lx52j;53p " 5lc5l53p ^^^ j + 27753x, (All) 

A-^ = (251 -52 + 53 + 251a; -522; + 253p- 53a; )53t + 277534, (A12) 



I 

where A is the corresponding RG flowing scale and 77 lous dimension is related to the self-energy feedback into 
is naturally the anomalous dimension defined in Eq. ([5|). the two- loop RG flow equations. 
As we have discussed previously in this paper, the anoma- 
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